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Abstract 


Collocation  implementation  of  the  Kleiser-Schumann’s  method  in  geometries  with  two  non¬ 
periodic  directions  is  shown  to  suffer  from  three  spurious  modes  -  line,  column  and  checkerboard 
-  contaminating  the  computed  pressure  field.  The  comer  spurious  modes  are  also  present  but  they 
do  not  affect  evaluation  of  pressure  related  quantities.  A  simple  methodology  in  the  inversion  of  the 
influence  matrix  will  efficiently  filter  out  these  spurious  modes. 


I .  This  research  was  partly  supported  by  the  National  Aeronautics  and  Space  Administration 
under  NASA  Contract  No.  NASl-19480  while  the  first  author  (SB)  was  in  residence  at  the 
Institute  for  Computer  Applications  in  Science  and  Engineering  (ICASE),  NASA  Langley  Re¬ 
search  Center,  Hampton,  VA  23665.  Authors  are  grateful  to  D.  Funaro  and  C.L.  Streett  for  their 
helpful  discussions. 


1.  Introduction 


Most  existing  numerical  solutions  of  incompressible  Navier-Stokes  equations  in  three  di¬ 
mensions  employ  the  primitive-variable  formulation.  Here,  the  velocity  field  and  pressure  can  not 
be  approximated  independently  and  must  satisfy  a  compatibility  condition.  Approximating  both  ve¬ 
locity  and  pressure  by  polynomials  of  the  same  degree  will  result  in  some  spurious  modes  for  the 
pressure  [  1  ].  In  a  spectral  collocation  implementation,  these  spurious  modes  can  be  characterized 
as  the  spurious  components  of  the  pressure  field  whose  discrete  gradient  at  the  interior  collocation 
points,  where  discretized  momentum  equations  are  satisfied,  is  zero.  Such  pressure  components 
have  no  effect  on  velocity  and  are  therefore  left  uncontrolled  by  the  discretized  governing  equations. 

Existence  of  spurious  pressure  modes  in  a  spectral  simulation  was  first  pointed  out  by  Mor- 
choisne  [2|.  In  coupled  spectral  implementations  where  the  continuity  equation  is  discretized  direct¬ 
ly,  theoretical  analysis  of  the  spurious  modes  is  on  a  firm  footing  [1].  When  a  non-staggered  grid 
is  employed  it  has  been  shown  that  a  fully  periodic  problem  has  no  spurious  modes,  discounting  the 
arbitrariness  of  the  mean  value  of  pressure  in  an  incompressible  flow.  In  a  flow  with  one  non-period¬ 
ic  direction  pressure  has  one  spurious  mode,  and  problems  with  two  non-periodic  directions  have 
seven  {l-line,  l-colunw,  \ -checkerboard  and  4-corner)  spurious  modes  [1,3].  These  results  have 
been  shown  to  apply  to  the  following  spectral  implementations:  Galerkin  Legendre,  Legendre  tau, 
Galerkin  Chebyshev,  Chebyshev  tau  and  Chebyshev  collocation  with  Gauss-Lobatto  points.  In  a 
collocation  implementation,  the  spurious  modes  can  be  avoided  with  an  appropriately  staggered 
mesh  for  velocity  and  pressure.  Montigny-Rannou  and  Morchoisne  [4]  have  recently  described  an 
algorithm  for  the  two  non-periodic  problem  on  a  half-staggered  grid  which  involves  only  one 
(checkerboard)  spurious  mode.  The  spurious  modes  can  be  completely  avoided  in  a  two  non-period¬ 
ic  problem  with  a  fully-staggered  mesh  [5,  6|. 

A  second  class  of  numerical  algorithms  avoids  direct  solution  of  the  continuity  equation  by 
solving  a  Pois.son  equation  for  pressure  1 1  [.The  time-split  implementation  [7]  is  the  simplest  of  them 
all  but  it  suffers  from  non-zero  boundary  divergence,  non-zero  slip  velocity  and  time-splitting  errors. 
Kleiser-Schumann’s  influence  matrix  method  can  be  used  to  enforce  zero  boundary  divergence  and 
no-.slip  condition  1S,9|.  A  partial  implementation  of  the  influence  matrix  is  simpler,  but  will  result 
in  small  but  non-zero  interior  divergence  [lOj.  A  full  implementation  of  the  Klei.ser-Schumann’s 
method  with  tau  (or  c(dlocation)  corrcclnm  enforces  interior  zero  divergence  as  well.  Here  we  iden- 
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tify  the  spurious  modes  in  the  collocation  implementation  of  the  Kleiser-Schumann  method  [8,  9, 
11]  with  tau  correction  for  a  problem  with  two  non-periodic  directions  to  be  the  line,  column,  check¬ 
erboard  and  comer  modes.  A  simple  correction  procedure  which  will  automatically  filter  these  spu¬ 
rious  modes  is  obtained.  This  correction  procedure  is  applied  in  the  simulation  of  a  turbulent  square 
duct  flow  [11]  and  is  found  to  be  very  effective  in  eliminating  spurious  modes. 

2.  Spurious  Modes  for  KJeiser-Schumann’s  Method 

The  full  implementation  of  the  Kleiser- Schumann’s  method  with  collocation  correction  can 
be  given  by  the  following  equations  and  boundary  conditions: 
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where  NL  is  the  nonlinear  term  in  the  Navier-Stokes  equation  and  Vh  is  the  velocity  boundary  condi¬ 
tion.  The  above  equations  and  boundary  conditions  are  in  their  discretized  form,  therefore  the  sym¬ 
bols  V,  V  •  and  repre.seni  discrete  gradient,  divergence  and  Laplacian  operators,  and  Q  and  dQ 
represent  interior  and  boundary  collocation  points.  vSince  the  momentum  equation  is  satisfied  only 
in  the  interior,  B  rcpre.sents  the  boundary  momentum  residual  and  is  nonzero  only  on  the  boundary. 
Combining  equation  (la)  and  the  di.screte  divergence  of  equation  (lb)  one  obtains 

— ^  provided  the  discrete  divergence  of  the  discrete  gradient  operator  is 

identically  equal  to  the  di.screte  Laplacian  operator.  The  above  equation  for  the  velocity  divergence 
along  with  the  boundary  condition  given  by  equation  Id,  results  in  a  divergence  free  velocity  field 
both  in  the  interior  and  on  the  boundary. 

By  definition,  each  spurious  mode  is  a  valid  .solution  to  the  di.scretized  governing  equations 
and  appropriate  boundary  conditions.  The  spurious  modes  have  a  non-zero  contribution  to  pre.ssure 
but  have  no  effect  on  velocity,  therefore  V,,,  =  0,  where  subscript  “sp”  .stands  for  the  .spurious 


mode.  The  spurious  pressure  components  therefore  satisfy 


VVsp  =  ^  '  Bsp  in 

0  =  -  Vpsp  +  Bsp  in  Q,  di2 


Bsp  =  0  in 

where  Bsp  is  the  corresponding  spurious  boundary  momentum  residual.  Any  non-trivial  solution  to 
the  above  linear  equations  represents  a  spurious  mode,  which  when  added  to  the  true  solution  will 
still  satisfy  the  discretized  Navier-Stokes  equations  and  boundary  conditions. 

Eight  solutions  to  the  above  equations  can  be  identified.  First  of  which  is  the  non-spurious 
solution,  =  constant  and  Bsp  =  0,  which  indicates  that  pressure  is  evaluated  only  up  to  an  arbitrary 
additive  constant  in  incompressible  flows.  The  first  two  spurious  modes  are  the  line  and  column 
modes 


Psp  -  T’/VjcW  .  Bi 


(±  lAiVx)2  at  jc  =  ±  1 
0  otherwise 


Pisp  -  B^sp  ~  ® 


Bsp  ~  ^A/yCv)  >  B- 


i±  lAyVy)^  aty  =  ±  1 
0  otherwise 


Bup  =  B,sp  =  0 


Here  iNx+ 1 )  and  (A/y+l )  are  the  number  of  points  along  the  non-periodic  Chebyshev  directions  and 
the  third  direction  is  at  most  periodic.  Bj ,  B2  and  Bj  are  th'’  three  components  of  the  boundary  mo¬ 
mentum  residual  and  only  the  normal  component  is  nonzero.  Tf^x  and  T^y  are  Chebyshev  polyno¬ 
mials  of  the  highest  degree  along  x  and  y.  The  third  spurious  mode  is  the  checkerboard  mode,  with 


Bsp  — 


=1: 


±  1)"'^(N.x)2  at  ;c  =  ±  1 


otherwise 


(±  lAAy)2  aty=±l 

0  otherwise 


B^sp  =  0 


These  three  spurious  modes  have  no  variation  in  the  periodic  z  direction  and  therefore  contaminate 
only  the  zero**’  mode  along  the  z  direction.  The  other  four  spurious  modes  are  the  comer  modes  and 
each  of  them  can  have  arbitrary  variation  along  the  z  direction.  For  example,  let  /j  j(z)  be  the  arbi¬ 
trary  variation  along  the  jr=l,  y=l  corner.  Then  the  comer  mode  corresponding  to  this  comer  can 
now  be  written  as 
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p  _  f/i,i^z>  at  Jf  =  V  =  1  _  f/i.i{2)  Df^^ix)  at  y  =  +1 

'  5p  I  *  ^  \SP  ^  I 

[0  otherwise  [0  otherwise 

B  ^  at  jc  =  +  I  _  _  at  x  =  v  =  1 

0  otherwise  ^  otherwise 

where  Df^x(x}  is  the  discrete  derivative  of  the  polynomial  which  collocates  to  zero  at  all  points  except 
at  ^1 .  The  comer  modes  for  the  other  three  comers  can  be  written  similarly.  The  comer  spurious 
modes  simply  reflect  the  fact  that  in  a  collocation  implementation  the  pressure  along  the  four  corner¬ 
lines  never  enter  into  the  computation  and  therefore  their  values  remain  unspecified.  This  arbitrari¬ 
ness  associated  with  the  corner  pressure  values  is  relatively  innocuous,  since  they  are  not  required 
in  computing  pressure  related  quantities  of  interest.  On  the  other  hand,  the  other  three  spurious 
modes  are  buried  in  the  numerically  computed  pressure  and  need  to  be  filtered  out. 

3.  Filtering  Procedure 


Implementation  of  the  Kleiser-Schumann’s  method  involves  the  construction  of  an  influence 
matrix.  The  solution  of  the  pressure  Poisson  equation  (equation  la)  requires  the  knowledge  of  pres¬ 
sure  boundary  conditions  {ph)  and  boundary  momentum  residuals  {Bh)  at  the  {2Nx+2Ny-A)  points, 
excluding  the  corner  points.  Note  that  in  the  three,  tine,  column  and  checkerboard  modes  only  the 
normal  component  of  the  boundary  momentum  residual  is  nonzero.  In  essence,  a  total  of 
(4Nx+4Ny~?>)  unknown  boundary  pressure  and  normal  momentum  residuals  are  required  in  evaluat¬ 
ing  the  pressure  field.  These  unknowns  are  evaluated  by  requiring  that  continuity  and  normal  mo¬ 
mentum  (with  the  residual)  are  satisfied  at  the  boundary  points  excluding  the  corner  points.  This  pro¬ 
vides  the  necessary  (4Nx+4Ny-^)  linear  equations  for  the  unknown  quantities.  These  equations  can 
be  cast  into  the  following  matrix  form,  Ax  =  R,  where  ^  is  the  influence  matrix,  j[  is  the  unknown 

vector  of  boundary  pressure  and  normal  momentum  re.sidual  and  R  is  the  right  hand  side.  In  a  three- 
dimen.sional  problem,  a  Fourier  tran.sform  along  the  periodic  z  direction  will  result  in  one  influence 
matrix  for  each  Fourier  wavenumber  k^.  Invcrtibility  of  the  influence  matrix  is  closely  related  to  the 
pre.sence  of  .spurious  modes.  In  particular  the  influence  matrix  corresponding  to  k^=i)  suffers  from 
the  constant,  line,  column  and  checkerboard  mode.s.  This  influence  matrix  therefore  has  four  zero 
eigenvalues.  In  the  numerical  computation  of  a  turbulent  square  duct  flow  it  was  observed  that  only 
for  A-=0  were  four  eigenvalues  of  the  influence  matrix  of  the  order  10“’“,  while  the  lowest  of  all  other 
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eigenvalues  were  around  10'^.  The  eigenvectors  corresponding  to  these  eigenvalues  are  the  corre¬ 
sponding  spurious  boundary  pressure  and  normal  component  of  the  momentum  residual. 


As  suggested  by  Tuckerman  [9],  the  non-invertibility  of  the  influence  matrix  can  be  easily 
overcome  by  constructing  a  related  matrix  ^  where  ^  and  *  are  the  eigenvec¬ 

tor  matrix  of  the  original  influence  matrix  and  its  inverse,  respectively,  and  A’  is  a  diagonal  matrix 
with  the  eigenvalues  along  the  diagonal  and  with  the  zero  eigenvalues  replaced  by  some  non-zero 

constant.  The  influence  matrix  now  becomes  invertible,  i.e.,  ==  M  1/A'  and  the  result¬ 


ing  pressure  and  boundary  momentum  residuals  yield  a  divergence-free  flow  field  independent  of 
the  constant  that  replaces  the  zero  eigenvalue.  Let  the  eigenvalue  of  the  original  influence  matrix 
be  zero  and  be  replaced  by  a  constant  Cp.  Let  the  corresponding  p‘^  eigenvector  be  Mip,  which  is  a 
vector  of  boundary  pressures  and  momentum  residuals  corresponding  to  a  linear  combination  of  the 


spurious  modes.  Contribution  of  this  p'^  mode  to  the  unknown  vector  x,  is  then 


,  where  bp  is 


the  projection  of  the  right  hand  side  along  the  eigenvectors,  given  by  M~  ‘  Rj.  One  simple  way  to 
filter  the  four  spurious  modes  will  then  be  to  set  the  arbitrary  constant  Cp  to  be  infinity.  In  other  words, 
in  the  evaluation  of  ^  ',  one  over  the  zero  eigenvalue  is  replaced  by  zero.  This  filtering  procedure 

was  implemented  in  the  computation  of  turbulent  flow  in  a  square  duct  [11]  and  it  essentially  re¬ 
moved  all  high  frequency  checkerboard  type  oscillations.  It  should  be  pointed  out  that  such  oscilla¬ 
tions  can  be  observed  in  the  k.^0  modes  as  well,  especially  at  high  Reynolds  numbers  when  resolu¬ 
tion  is  only  marginal.  These  oscillations  are  not  due  to  the  spurious  modes,  but  are  controlled  energy 
in  the  high  wavenumber  modes. 


4.  Spurious  Modes  in  Partial  and  Time-Split  Methods 


Spurious  modes  for  the  partial  implementation  of  the  Kleiser-Schumann’s  method  without 
the  collocation  correction  can  be  analyzed  in  similar  fashion.  The  spurious  modes  satisfy  the  follow¬ 
ing  equations. 
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^^Psp  =  0 


in  Q 


0  =  -  Vpsp  +  Bsp 


Bsp  —  0 


in  Q,dQ 


in  Q 


(6a-6c) 


It  can  be  easily  seen  that  the  only  admissible  solution  for  the  above  set  of  equations  is  the  non-spuri- 
ous  constant  mode  and  there  are  no  spurious  modes  present.  This  is  confirmed  in  the  numerical  simu¬ 
lation  by  observing  that  the  influence  matrix  for  the  mode  has  only  one  zero  eigenvalue  and 
the  corresponding  constant  mean  pressure  can  be  set  to  zero  by  replacing  the  zero  eigenvalue  by  in¬ 
finity. 

Spurious  modes  in  the  time  split  implementation  will  depend  on  the  exact  boundary  condi¬ 
tions  employed  for  the  intermediate  star-level  velocities  (V*)  and  the  pressure  Poisson  equation.  Fol¬ 
lowing  Streett  and  Hussaini  [7],  if  we  employ  [(2Vp(r)  —  Vp{t  —  dt))  •  r]  as  the  boundary  condition 

for  the  tangential  components  of  the  intermediate  star-level  velocity,  zero  penetration  for  the  normal 
velocity  component  and  zero  Neumann  boundary  condition  for  the  pressure,  then  we  have  the  fol¬ 
lowing  equations  satisfied  by  the  spurious  modes: 


^sp  ~  Vpsp 
/it  ^P 

^  sp  '  V  ~  ^Psp  '  tj  —  0 
V*sp’Z  =  {l^Pspit)  -  Vpspit  -  dt))  ■  T 


in  Q 


in  Q 


in  dQ 


in  dQ 


(7a-7d) 


where  ij  and  t  are  direction  normal  and  tangential  to  the  boundary.  The  analysis  of  the  spurious 
modes  is  more  complicated  and  also  depends  on  the  initial  tangential  pressure  gradients  on  the 
boundary.  With  careful  choice  of  initial  conditions,  the  no  penetration  and  pure  Neumann  boundary 
conditions  will  guarantee  no  spurious  components. 

5.  Conclusion 

Collocation  implementation  of  the  Kleiser-Schumann’s  method  in  geometries  with  two  non¬ 
periodic  directions  have  three  spurious  modes -/me,  column  and  checkerboard -contaminating  the 
computed  pressure  field.  The  corner  spurious  modes  are  also  present  but  they  do  not  affect  evalua¬ 
tion  of  pressure  related  quantities.  The  three  spurious  modes  can  be  easily  filtered  out  by  replacing 
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the  zero  eigenvalues  of  the  influence  matrix  with  infinity  before  solving  for  the  unknown  boundary 
pressure  and  momentum  residuals.  Partial  implementation  of  the  Kleiser-schumann’s  method  with¬ 
out  collocation  correction  admits  no  spurious  modes.  Spurious  modes  can  also  be  avoided  in  time- 
split  implementations. 
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